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_ I  To  ^compute  the  conjugate  depth  at  a  hydraulic  jump  in  a  circular  open  channel.. 

References :  (l)  Chow,  Ven  Te,  Open  Channel  Hydraulics,  McGraw-Hill,  1959* 

pp.  53-56. 

(2)  Schied,  Francis,  Numerical  Analysis,  Schaum's  Outline  Series, 
O  McGraw-Hill,  1968,  pp  315-317. 
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^>The  program  is  written  in  G635  time-share  series,  FORTRAN  IV,  and  is  part  of  a 
Conversationally  Oriented  Real-Time  Program-Generating  System  (CORPS).  The 
program  consists  of  a  main  program  and  three  subroutines.  The  main  program 
handles  all  I/O  requirements  with  the  subroutines  handling  the  computation.  -v 


0.  RROORAM  SPECIFICATIONS 


SEE  FOLLOWING  PAGE 
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D.  IOUIMINT  OKTAIU 


The  program  was  developed  and  is  operational  on  the  WES  G635,  Vicksburg,  MS. 
It  is  also  operational  on  HIS  66/80,  Macon,  GA,  and  Boeing  CDC,  Seattle,  WA. 


|  B.  MPUf-WTWT - - 

The  required  inputs  which  are  entered  at  execute  time  are:  discharge,  cfs; 
diameter,  ft;  number  of  flow  depths,  and  flow  depth,  ft.  The  output  consists 
of  the  given  data,  plus  the  flow  depth  and  corresponding  conjugate  depth  in 
tabular  form. 


F.  ADDITIONAL  ft  (MARKS 


Complete  documentation  of  this  program  is  available  from  the  Engineer 
Computer  Programs  Library,  Technical  Information  Center,  WES. 
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B.  PROGRAM  SPECIFICATIONS: 

Language:  ANSI  FORTRAN  (FORTRAN  IV) 

Solution  Requirements:  The  run  command 

RUN  WESLIB/ C0RPS/H6122 ,  R 

plus  the  input  variables  defined  in  (E). 

Method  of  Analysis:  Solves  the  momentum  equation  for  the  conjugate 
depth  by  application  of  Newton’s  Tangent  Method. 

Core  Requirements  G635:  11  K  words 

External  Storage:  None 

Restrictions:  Velocity  distribution  is  uniform  across  flow  section 
and  the  invert  slope  is  small  (i.e.,  <  10  degrees). 

General  Equation: 

A.y.  +  =  A_y_  +  ^ 


fl  gA1  V2  gA2 


where:  Q  is  discharge  (cfs),  g  is  acceleration  of  gravity 
2 

(32.2  ft/sec  ),  A^  and  Ag  are  the  respective  cross-section  areas 

2  -  — 

(ft  )  for  flow  depth  and  conjugate  depth,  and  y^  and  yg  are  the 

distances  (ft)  of  the  centroid  of  the  respective  water  areas  A^ 

and  Ag  below  the  water  surface. 

Range  of  Quantities:  Unlimited  for  practical  application. 


computed  to  ±0.001  ft  in  H6122. 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  Computer 
Program  Library  Standards  and  Documentation,  Appendix  B 

PART  I:  ENGINEERING  DESCRIPTION 

1.  PROGRAM  NUMBER:  722-F3-R0-6AH 

2.  TITLE:  H6122  -  Conjugate  Depth  in  a  Circular  Open  Channel 

3.  REVISION  LOG:  N/A 

U.  PURPOSE  OF  PROGRAM:  To  compute  conjugate  depth  at  a  hydraulic  Jump 
in  a  circular  open  channel. 

References: 


a.  Chow,  Ven  Te,  Open  Channel  Hydraulics.  McGraw-Hill,  195 9* 

PP  53-56. 

b.  Schied,  Francis,  Numerical  Analysis.  Schaum's  Outline  Series, 
McGraw-Hill,  1968,  pp  315-317. 


5.  STEP  SOLUTION: 


a.  The  inputs,  discharge  (Q)  cfs,  diameter  (D)  ft,  and  flow  depth 
(y1)  ft,  are  entered. 

b.  Computational  steps: 


(1) 

(2) 


Specific  force 

Specific  force 
flow  and  F1 


(F^)  at  y^  using  subroutine  H6125 
difference  (FD)  between  that  at  full  conduit 


FD 


irD" 


girD' 


2  +  T 


-  F, 


(3)  If  FD  <  0,  the  conjugate  depth  (y2)  >  D.  The  flow  depth 

(y^)  is  printed  with  the  message  CON  DEPTH  >  DIAMETER. 

Control  is  then  returned  to  the  main  program.  If  FD  >.  0, 
then  step  (k)  is  entered. 


(4)  Critical  depth  (yc)  using  subroutine  H6lUl 
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(5)  Determine  trial  conjugate  depth 
Ye 

y+  “  *- 


(yt) 


If  y 

*  n 


then 


If  yc  =  yi  ,  then  conjugate  depth  (y2)  is  equal  to  yc  ; 
yp  is  returned  to  the  main  program.  If  y  >  y.  »  then 
D  +  yc 


.  _  £ 

SPECIFIC  FORCE  CURVE 


(6)  Specific  force  (F  )  and  area  (A.  )  at  y.  using  sub¬ 
routine  H6125 


(7)  Define  function  F(y.)  to  be  used  for  Newton's  Tangent 
Method 


F(yt)  -  Ft 


F, 


(8)  First  derivative  of  F(y^)  with  respect  to  y^ 


At  “  ^2  *ytD  “ 


gA- 


1/2 

) 
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(9)  Conjugate  depth  (y2)  is  calculated  using  Nevton's  Tangent 
Method 

F(yJ 

y2  *  yt  "  ”d - 

Z  —  F(y  ) 

dyt  FKyt> 


(10)  If  y„  <  0.0  ,  then  y ■  *  —  ,  where  m  *  2. 3 .depending 

i  —  t  gin 

on  the  number  of  times  the  condition  yg  <  0.0  occurs. 

Transfer  is  to  step  (6)  and  the  process  is  repeated  until 
yg  >  0.0  .  Step  (11)  is  then  entered. 


(11)  If  y0  >  D  ,  then  y *  D  - 


D  -  yc 


,  where  m  *  2,3*. . . , 


depending  on  the  number  of  times  the  condition  y2  >  D  , 
occurs.  Transfer  is  to  step  (6)  and  the  process  is 
repeated  until  yg  <_  D  .  Step  (12)  is  then  entered. 

(12)  If  |yg  -  y^|  >.  .001  ,  then  y^  is  set  equal  to  yg  and  the 

procedure  transfers  to  step  (6)  and  is  repeated  until  the 
|y2  -  y^|  <  .001  .  When  |y2  -  yt|  <  .001  ,  then  the 

value  of  the  conjugate  depth  (y2)  is  returned  to  the 
main  program. 

c.  The  given  data,  plus  the  conjugate  depth,  are  printed. 

6.  ACCURACY ;  Governed  by  accuracy  of  input  data;  conjugate  depth  is 
computed  to  +.001  ft. 

7*  REMARKS:  Velocity  distribution  is  uniform  across  flow  section  and 
the  invert  slope  is  small  (i.e.,  <  10  degrees). 
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PART  II:  COMPUTER  FUNCTIONAL  DESCRIPTION 
REVISION  LOO:  N/A 
FUNCTIONAL  FLOW  CHART: 
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~[g6l22  |  | 

Jl 

(Specific  force  for  flov  depth  using  subroutine 
H&125;  specific  force  difference  betveen  force 
Ut  flov  full  end  flov  depth 


PRUT:  flov  depth 
end  Message  COM 
DEPTH  >  DIAMETER 


^RETURN  1^ 


Use  subroutine 
H6lUl  to  compute 
critical  depth 


Trial  conjugate  depth 
_<  critical  depth 


Trial  conjugate  depth 
>  critical  depth 


Starting  vith  trial  conjugate  depth  and  using 
Mevton'1  s  Tangent  Method  as  a  depth  correction. 

compute  conjugate  depth 

H6122 


3-  EQUIPMENT  AMP  OPERATING  SYSTEM:  The  program  was  developed  on  a 
G635  tlmeshare  system  in  which  input/output  equipment  consisted  of  a 
Model  33  remote  teletype.  It  is  now  operational  on  the  WES  G635* 
Vicksburg,  MS;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 

U.  INPUT  REQUIREMENTS:  The  required  inputs  are  entered  via  the  user's 
time-share  terminal  device  in  free  field  format.  All  input  cues  and 
reads  are  performed  in  the  main  program.  The  subroutines  handle  the 
computations.  Since  computations  are  done  in  the  subroutine,  and  sub¬ 
routines  H6125  and  H6lUl  are  subordinate  to  H6122,  the  necessary  inputs 
to  subroutine  H6122  are  passed  via  the  CALL  statement.  The  calling 
sequence  is: 

CALL  H6l22( arg1 , . . . ,arg^ ,  $N ) 

where : 

arg^^  -  discharge,  cfs 
argg  -  diameter,  ft 
arg^  -  flow  depth,  ft 
argjj  -  conjugate  depth,  ft 

$N  -  N  is  a  statement  number  for  the  nonstandard  return  due  to 
conjugate  depth  >  diameter 

All  arguments  cure  floating  point,  except  N,  which  is  integer. 

Arguments  1-3  are  inputs  and  arg^  is  output.  N  is  a  STATEMENT  NUMBER 

which  directs  the  return  from  the  subroutine  to  that  STATEMENT  NUMBER 

in  the  main  program. 
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5.  SECONDARY  STORAGE  INPUT:  None 

6.  INPUT  DATA  DESCRIPTION;  The  following  names  are  used  for  the  input 
variables  in  program  H6122. 

DISCH  -  discharge,  cfs  ( arg^ ) 

DIAMTR  -  diameter,  ft  (argg) 

N  -  number  of  flow  depths,  integer 

Y  -  flow  depth,  ft  (arg^) 

7.  OUTPUT  DATA  DESCRIPTION:  The  following  names  are  used  for  the 
output  variables  in  program  H6122. 

DEPTH2  -  conjugate  depth,  ft  (arg^) 

8.  PROGRAM  ERROR  MESSAGES: 

a.  In  the  main  program,  if  the  flow  depth  is  £  0  or  >  diameter  the 
following  message  is  printed: 

FLOW  DEPTH  6.00  ft  <  OR  =  0  OR  >  DIAMETER 

The  message  is  repeated  until  a  valid  flow  depth  is  entered. 

b.  In  subroutine  H6122,  the  value  of  the  flow  depth  is  printed 
and  then  the  message  CON  DEPTH  >  DIAMETER  is  printed  on  the 
same  line. 

For  example:  1.20  CON  DEPTH  >  DIAMETER 

Return  is  then  to  the  statement  number  in  the  main  program 
specified  by  N  in  the  CALL  statement. 

9.  VARIABLE  DEFINITIONS: 

N  -  number  of  flow  depths,  0  <  N  <_  25 

y  -  flow  depth,  ft;  dimensioned  max  25 

DEPTH1  -  flow  depth,  ft;  equal  to  Y(I)  for  I  =  1,...,N 
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CRIT  -  critical  depth,  ft;  returned  from  subroutine  H6l4l 

DEPTH 2  -  conjugate  depth,  ft 

DEPTHT  -  trial  conjugate  depth,  ft 

DIAMFR  -  diameter,  ft 

DISCH  -  discharge,  cfs 

F  -  function  of  trial  conjugate  depth  used  for  Newton1 s 

Tangent  Method,  ft3;  the  value  of  the  function  is  a 
specific  force  difference  between  that  at  trial  conjugate 
depth  and  flow  depth 

FP  -  first  derivative  of  function  (F)  with  respect  to  trial 

conjugate  depth,  ft2 

FORCE  -  specific  force,  ft^;  returned  from  subroutine  H6125 

FI  -  working  storage  to  hold  FORCE  at  a  constant  when  FORCE  is 

returned  from  H6125  as  the  specific  force  at  flow  depth, 
ft  3 

FMAX  -  specific  force  difference  between  that  at  full  conduit 
flow  and  flow  depth,  ft 3 

p 

G  -  acceleration  of  gravity,  32.2  ft/sec 

HFILB  -  five  character  name  of  program;  passed  to  WES LIB  count 
routine  HACCT 

LQZ  -  equal  1,  execute  all  input  cues  and  reads;  equal  2,  call 
WESL1B  routine  RERUN  and  enter  only  desired  inputs 

LQX  -  equal  1,  print  instructions  from  RERUN;  equal  3,  no  print 

JKL  -  direct  return  from  RERUN  to  desired  input  read 

KXK  -  total  number  of  inputs;  passed  to  RERUN 

M  -  powers  of  2  to  converge  Newton's  Tangent  Method  if  trial 

conjugate  depth  £  0  or  >  diameter 

PI  -  constant  for  ir;  value  *  3.1U159265 
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WK1 


-  working  storage  =  T°P  from  subroutine 

H6125 


ZZZZZ  -  character;  equal  RE,  rerun;  equal  ST,  stop 
10.  EXAMPLE  CASE:  Compute  the  conjugate  depth  for  6  given  flow  depths. 


a.  Input  data: 

Discharge  (DISCH)  *  100.00  cfs 
Diameter  (DIAMTR)  =  5-00  ft 
Number  of  flow  depths  (N)  =  6 

Flow  depths  (Y(l),  for  I  =  1,N)  =  1.00,  1.43,  2.50,  4.00,  4.50, 
and  5-00  ft 

b.  Output : 


INPUT  H6 122-CON JUGATE  DEPTH(S)  IN  A  CIRCULAR  OPEN  CHANNEL 

AA-ENTER  DISCHARGE, CFS . 

=  100 

AB-ENTER  DIAMETER. FT. 

=5 

AC-ENTER  THE  NUMBER  OF  DEPTH(S)  FOR  WHICH  THE  CONJUGATE  DEPTH  IS  TO  BE 
CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS. 

=6 

AD-ENTER  THE  6  DEPTH(S)  SEPARATED  BY  COMMAS. 

=1,1.43,2.5,4,4.5,5 


OUTPUT  H6 122-CONJUGATE  DEPTH(S)  IN  A  CIRCULAR  OPEN  CHANNEL 

DISCHARGE  =  100.00  CFS 

DIAMETER  =  5.00  FT 


FLOW  CONJUGATE 
DEPTH  DEPTH 

(FT)  (FT) 

1.00  CON  DEPTH  >  DIAMETER 

1.43  CON  DEPTH  >  DIAMETER 

2.50  3.23 

4.00  1.97 

4.50  1.72 

5.00  1.52 


ENTER  RERUN  OR  STOP 
=ST0P 
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REF:  ER  1110-1-10  -  ENGINEERING  AND  DESIGN  -  Engineering  and  Computer 
Program  Library  Standards  and  Documentation,  Appendix  C 

PART  Ills  FILE  DOCUMENTATION 

1.  REVISION  LOG:  N/A 

2.  TITLE:  H6122  -  Conjugate  Depth  in  a  Circular  Open  Channel 

3.  SOURCE  PROGRAM  LISTINGS:  See  pages  12-15 

4.  NUMERICAL  AND  LOGICAL  ANALYSIS:  Solves  the  momentum  equation  for 
the  conjugate  depth  by  application  of  Newton's  Tangent  Method. 

5.  SUBROUTINES  NOT  DOCUMENTED  IN  ABSTRACT:  The  following  subroutines 
are  used  in  program  H6122. 

a.  SUBROUTINE  H6125  (DISCH,  DIAMTR,  DEPTH,  A,  WK1,  FORCE) 

The  subroutine  statement  for  H6125  as  documented  in  program  H6125  is: 

SUBROUTINE  H6125  (DISCH,  DIAMTR,  DEPTH,  ENERGY,  FORCE) 

The  argument  list  was  changed  to  (a.)  to  facilitate  the  solution  of 
program  H6122.  This  change  has  no  effect  on  H6125  as  documented. 

b.  SUBROUTINE  H6lUl  (DISCH,  DIAMTR,  DEPTH,  CRIVEL) 

Complete  documentation  of  these  subroutines  is  available  from  the 
Engineer  Computer  Program  Library,  Technical  Information  Center,  WES. 

6.  Miscellaneous :  The  program  is  part  of  the  CORPS  computer  system. 
CORPS  is  an  acronym  standing  for  Conversationally  Oriented  Real-Time 
Program-Generating  System.  The  program  is  now  operational  on  the  WES 
G635,  Vicksburg,  KB;  HIS  66/80,  Macon,  GA;  and  Boeing  CDC,  Seattle,  WA. 
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The  source  listing  on  page  12  contains  the  first  line  run  command  and 
brief  for  H6122.  This  first  line  run  command  runs  the  binary  H6122B 
of  the  source  listing  on  pages  13-15  (Fortran  source  of  H6122)  and 
attaches  the  WESLIB  routines  HACCT  and  RERUN. 
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0 0 0 X « ORUN  klESL IB/C0RPS/H6 122B , R ; UESL IB/RERUN , R ; MESL IB/HACCT . R 

0800  63THIS  PROGRAM  COMPUTES  CONJUGATE  DEPTH(S)  IN  A  CIRCULAR  CHANNEL. 

0805  36VELOCITY  DISTRUBUTION  COEF.  =  UNITY. 

0810  61  INPUTS  REQUIRED  ARE  DISCHARGE-CFS.  DIAMETER-FT,  AND  FLOW 
0815  09DEPTH-FT . 

0820  53  OUTPUT  INCLUDES  THE  GIVEN  DATA  AND  CONJUGATE  DEPTH 
0825  58  IF  CONJUGATE  DEPTH  >  DIAMETER. THEN  THE  VALUE  OF  THE  FLOW 
0830  55DEPTH  IS  PRINTED  AND  THE  MESSAGE  CON  DEPTH  >  DIAMETER. 
0999X06FINISH 
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00001X0RUN  *=»/C0RPS/H6122B(N0G0) 

10000  CHARACTERM5  HFILE 
10010  DIMENSION  Y(25) 

10020  HFILE=5HH6122 
10030  LQZ*1 ; LQX=1 
10040  15000  CONTINUE 
10045  PRINT  10050 

10050  10050  FORMAT (/"INPUT  H6 1 22-CONJUGATE  DEPTH(S)  IN  A  CIRCULAR  OPEN 
100558CHANNEL"//) 

10060  CALL  HACCT (HFILE) 

10070  GO  TOU5003, 15016), LOZ 
10000  15003  PRINT  15004 

10090  15004  FORMAT ("AA-ENTER  DISCHARGE. CFS."> 

10100  15005  READ.DISCH 

10110  GO  T0(15007, 15016), LOZ 

10120  15007  PRINT, "AB-ENTER  DIAMETER, FT 

10130  15008  READ.DIAMTR 

10140  GO  T0(15010, 15016), LOZ 

10150  15010  PRINT  15111 

10155  15111  FORMAT("AC-ENTER  THE  NUMBER  OF  DEPTH(S)  FOR  UHICH  THE  CONJUO 
101578ATE  DEPTH  IS  TO  BE"/"CALCULATED.  MUST  NOT  EXCEED  25  DEPTHS.") 

10160  15011  READ.N 

10170  GO  T0(15013, 15016). LOZ 

10180  15013  PRINT  15113. N 

10185  15113  FORMAT ( "AD-ENTER  THE  ",I2,"  DEPTH! S)  SEPARATED  BY  COMMAS.") 
10190  15014  READ, (Y(I>. 1=1, N) 

10200  GO  T0(15019. 15016), LOZ 

10210  15016  KKK=4 

10220  CALL  RERUN(KKK.LQX.JKl) 

10230  GO  T005005, 15008, 15011, 15014, 15019), JKL 
10240  15019  PRINT  15020, DlSCH.DIAMTR 

10250  15020  FORMAT(/"OUTPUT  H6 122-CONJUGATE  DEPTH(S)  IN  A  CIRCULAR  OPEN 
102558  CHANNEL "//"DISCHARGE  =  ",F13.2,"  CFS"/"DIAMETER  =  ",F13.2,"  FT"// 
1026084X. "FLOW", 5X,"CONJUGATE"/4X. "DEPTH", 6X."DEPTH"/4X, "(FT )".7X, "(FT)" 
102708) 

10280  DO  15027  1*1, N 
10290  DEPTH1=Y( I ) 

10300  CALL  H6 1 22 (DISCH,DIAMTR,DEPTH1,DEPTH2, 015027) 

10310  PRINT  15028. DEPTH1 .DEPTH2 

10320  15027  CONTINUE 

10330  15028  FORMAT (F8.2.F11.2) 

10340  PRINT,"  " 

10350  LQZ*2 

10360  CHARACTER  ZZZZZZ«2 

10370  16000  PRINT,  "ENTER  RERUN  OR  STOP" 

10380  READ  16001,  ZZZZZZ 
10390  16001  FORMAT (A2) 

10400  IF (ZZZZZZ . EO . 2HRE)  60  TO  15000 
10410  I F ( ZZZZZZ . EG . 2HST )  GO  TO  20000 
10420  PRINT, "ERROR  *«*  RETYPE" 
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10430  GO  TO  16000 
10440  20000  STOP; END 

20000  SUBROUTINE  H6122( DISCH * DIAMTR* DEPTH1 *  DEPTHS* k ) 

20010  G=32.2*PI=3. 141 59265 ;N=2 

20020  CALL  H6 1 25  (  DI  SCH  *  DIAMTR  *  DEPTH  1  *  AREA  *I4K1 » FORCE) 

20030  Fl=FORCE 

20040  FMAX=PlxDIAMTRK»3/6.+4.XDI$CHX«2/(GKPI»DIAMTRXK2)-Fl 
20050  IF(FMAX)  10050 *  1 0060 . 10080 
20060  10050  IF(DIAMTR-DEPTH1 )  ,10060* 

20070  PRINT  10060 » DEPTH 1  * RETURN  1 

20060  10060  FORMAT (F6. 2 *"  CON  DEPTH  >  DIAMETER") 

20090  10080  CALL  H6141 ( DISCH * DI AMTR * CRIT *CV) 

20100  I F< CRIT -DEPTH 1 )  10100*10110*10120 

20110  10100  DEPTHT-CRIT/N  *N=2#N ;GO  TO  10130 

20120  10110  DEPTH2=CRIT ;RETURN 

20130  10120  DEPTHT=DIAMTR-(DIAMTR-CRIT)/N;N=2XN 

20140  10130  CALL  H6125C DISCH * DIAMTR * DEPTHT* AREA *WK1 *  FORCE) 

20150  F=F0RCE-F1 

20160  FPsAREA-2.*DISCH#x2*WKl/(G*AREA*X2) 

20170  DEPTH2 -DEPTH T-F/FP 

20180  IF<  DEPTHS . LE . 0 . )  GO  TO  10100 

20190  IF(DEPTH2.GT. DIAMTR)  GO  TO  10120 

20200  IF( AB$<  DEPTH2-DEPTHT )  •  LT . . 001 )  RETURN 

20210  DEPTHT -DEPTH2*G0  TO  10130 

20220  END 

30000  SUBROUTINE  H6125< DISCH * DIAMTR *  DEPTH , A ,WK1 *  FORCE) 

30010  G=32.2;WK1=5QRT<DEPTH*DIAMTR-DEPTH**2) 

30020  WK2=2.*DEPTH-DIAMTR;WK3-ARC0S<-WK2/DIAMTR) 

30030  A-WK1XWK2/2 . +WK 3# DIAMTR* #2/4 • 

30040  AY=UKl#(DEPTH##2/3.-DEPTHKDIAMTR/3.*DIAMTR#»2/4.)4DIAMTR#*2#WKS«U 
300501K3/8. 

30060  ENERGY=DEPTH  ♦DI$CH#»2/(2.#G#AK#S) 

30070  F0RCE=AY*DISCH##2/(G*A) 

30080  RETURN 
30090  END 

40000  SUBROUTINE  H6141 (DISCH* DIAMTR, DEPTH*CRIVEL ) 

40010  DEPTH=DIAMTR/2. 

40020  6-32 . 2 

40030  SFAC=DISCH/SQRT(G> 

40040  5  CALL  H6 1 25 ( DISCH » DIAMTR*  DEPTH* AREA* MK1 *  FORCE) 

40050  SECFAC=AREA#SQRT(AREA/(2.*UK1)) 

40060  DY=DEPTH/<  SFAC/SECFAC) -DEPTH 

40070  IFCDY.GT.-0.01. AND, DY. LT . 0 . 01 )  GO  TO  10 

40080  DY=DY/2. 

40090  J)EPTHC  =  DEPTH-DY 

40100  IFCDEPTHC. LE. 0.0. OR. DEPTHC.GE. DIAMTR)  DEPTHC*DEPTH/2 . 

40110  DEPTH=DEPTHC 
40120  GO  TO  5 

40130  10  CRIVEL=DISCH/AREA 
40140  RETURN 
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